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Abstract. The characteristic double-bumped spectral energy distribution (SED) of blazars is explained 
by either leptonic or hadronic models. In the former, Inverse Compton emission dominates the emission of 
the high energy bump, while proton synchrotron emission and proton-gamma interactions dominate it in the 
latter. 

We present a new stationary lepto-hadronic code that evaluates both the leptonic and the hadronic 
interactions. Apart from the modelling of the SED produced in a leptonic or hadronic model, the code 
permits the study of interesting mixed lepto-hadronic scenarios, where both processes contribute significantly 
to the high energy bump. 

A first application to data from the high frequency peaked BL Lac object PKS 2155-304 is discussed. 
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1 Introduction 



The spectral energy distribution (SED) of blazars is characterised by two non-thermal emission bumps, the first 
one peaking in the optical or in X-rays, the second one peaking in 7-rays. The origin of this emission is generally 
ascribed to a region of high density, filled with a magnetic field and moving in a relativistic jet aligned to the 
line-of-sight. 

While there is a general consensus that the low energy bump is due to synchrotron emission from electons 
and positrons in the emitting region, the origin of the high energy emission is interpreted differentely in leptonic 
or hadronic models. In the first scenario, the emission is entirely dominated by leptons, which produce the 
high energy bump through inverse Compton scattering off the synchrotron photons themselves (Synchrotron- 
Self-Compton, SSC) or off an external photon field. In the hadronic scenario the high energy bump is usually 
attributed entirely to proton synchrotron emission or photo-hadronic interactions. 

We have developped a lepto-hadronic code that reproduces both "extreme" scenarios and also allows the 
study of mixed scenarios with characteristic signatures in the high energy band. 



2 Code development 

2.1 Leptonic processes 

The framework of our model is the stationary one-zone SSC code developped by Katarzynski et al. (2001). A 



spherical region (characterised by its radius R) moves in the relativistic jet with Doppler factor S : it is filled 
with a homogenous magnetic field B and a primary electron population Tie{j—E/{mc'^)), which is parametrised 
by a broken power law (defined by the two slopes ai,2, the Lorentz factors 7e,mm , le,break , le,max and the 
normalization factor K^, at 7 = 1). 



The leptonic part of the original code has been improved as follows : 
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The synchrotron emission is evaluated computing the exact integration over pitch angles. The approxima- 



tion used by Katarzyriski et al. ( j2001j is valid only in the energy range of interest if the emitting particles 



are e 



The internal absorption due to 7-7 pair production is evaluated using the parametrization for the cross- 



section given by Aharonian et al. (20081, instead of a i5-function. The pair injection rate from 7-7 inter- 



actions is computed following Aharonian et al. (1983). A stationary distribution of secondary pairs 



(l*** generation), taking into account synchrotron cooHng, is evaluated following Inoue & Takahara (1996). 
Synchrotron emission from the stationary pair distribution is evaluated as well. 

• The absorption induced by the extra-galactic background light (EBL) on TeV photons is evaluated using 
the EBL model by Franceschini et al. (2008). 



2.2 Hadronic processes 

The spherical emitting region is filled with a primary proton population np{-f), parametrized by a power law 
with slope ap, normalization factor Kp at 7 = 1 and minimal and maximal Lorentz factor jp_min and jp^max- 
The proton synchrotron emission is corrected for internal 7 — 7 absorption and the associated first generation 
pair spectrum is evaluated as described above. 



In hadronic scenarios, an important role is played by photo-hadronic interactions. They have been computed 
using the Monte Carlo code SOPHIA ( |Mucke fc Protheroe||2001[ ). 

The target photon field is composed by the synchrotron radiation from primary electrons and from the pro- 
tons themselves. The energy of the interacting protons is corrected for synchrotron losses following Miicke & 



Protheroe (2001). The SOPHIA code is then called for 10 sampled proton energies. The distributions of the 
generated particles (7, e^,p, n, z^e./^, t'e./i) are summed and normalized to the number of protons in the blob 
suffering p-7 interactions. 

The spectra of the generated e* are corrected for radiative cooling to arrive at a steady-state solution. Syn- 
chrotron emission from secondary leptons and from the first generation pair spectrum from 7-7 interactions is 
then evaluated. 

/X* can emit synchrotron radiation before decaying into e*. They are retrieved from SOPHIA, and, in a first ap- 
proximation, their synchrotron emission is evaluated without considering the competition between synchrotron 
photon emission and decay into e^. 



3 Application to PKS 2155-304 data 

As a first test of our code, we have reproduced the SED of a bright high frequency peaked blazar for different 
scenarios. Fig. [T] shows an application to the data from the 2008 multiwavelength campaign on the blazar PKS 

). This data set represents at this time one of the most complete simultaneous 
SEDs of a TeV blazar: the shape of the high energy bump is fully characterised by the Fermi and H.E.S.S. 
data. PKS 2155-304 has been found to be in a low activity state during the campaign. In the interpretations 
presented here, the emission from the blob explains only the X-ray to 7-ray emission, assuming that from radio 
to visible light the emission from the extended jet dominates. 

We report the modelling of the SED in three different scenarios : a pure SSC leptonic model, a proton 
synchrotron dominated hadronic model and a third, mixed, scenario in which both the synchrotron emission 
from protons and the inverse Compton component contribute to the high energy bump. 

The modelling of the SED for the three different cases is shown in Fig. [l} and the model parameters are 
given m Table [1] The first scenario (top left in Fig. 1), is a standard SSC model, assuming a low magnetic field 
(75 mC), far from equipartition; the synchrotron emission from secondary particles from 7-7 pair production is 
more than three order of magnitude below the primary synchrotron emission. The second scenario (top right 
in Fig. 1) is a proton synchrotron dominated hadronic model, in which the inverse Compton component is 
almost three order of magnitude lower than the proton synchrotron emission; /x-cascade emission contributes 
(even though at ~ 1% only) to the high energy bump; the magnetic field value is much higher (80 G) than 
in the SSC case, in order to have a significant contribution from protons; the presence of protons significantly 



2155-304 (Aharonian etari|2009 
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increases the luminosity of the jet (equal to 10''^ ergs s^^) . The last case (bottom in Fig. 1) is an example of a 
mixed lepto-hadronic scenario, in which both inverse Compton and proton synchrotron emission contribute to 
the high energy bump; the magnetic field value is small (0.22 G), but the jet luminosity is high, given the high 
value of the emitting region size (ten times the value assumed for the proton synchrotron model). 




Fig. 1. Modelling of PKS2155-304 in three different scenarios. Colour code is as follow : primary electron synchrotron 
emission - red; inverse Compton emission - green; proton synchrotron emission - blue; muon synchrotron emission - violet. 
Dotted lines represent the synchrotron emission from the secondary pairs from 7-7 pair production. Top left: pure SSC 
model; top right: proton synchrotron dominated model; bottom : mixed lepto-hadronic scenario. 



4 Conclusions 

The code presented here permits to reproduce the observed blazar SED in different scenarios, spanning a wide 
parameter space: a leptonic or a hadronic solution can be found assuming different physical conditions (magnetic 
field, density and distribution of primary particles) in the blob. In addition, mixed lepto-hadronic scenarios 
naturally arise in this framework. 

The increasing data quality at high energies, provided by Fermi in the GeV range, and by current and 
planned Cherenkov telescopes in the TeV range, will hopefully help to evaluate contributions from both leptons 
and hadrons. The next step towards a more realistic emission model would be a time-dependent lepto-hadronic 
code, which could take into account constraints from the observed variability from blazars. 
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SSC Proton synch Mixed 



9 1° 1° 1° 

6 30 30 30 

7e,min 5 X 10^ 1 X 10^ 1 X 10^ 

-fe,break 1-2 X 10^ 4 X 10^ 6 X 10^ 

le,max 1 X 10^^ 5 X 10^ 5 X 10^^ 

ae,i 2.4 2.2 2.0 



ae,2 

Up 



4.32 4.5 4.25 



Ke 7 X lO'' 3.1 X 10^ 6.7 x 10^ 

5.4 X 10-3 4.8 X lO-'^ 1.3 x 10"^ 
7p,mzn - 1 X 10^ 1 X IQ^ 

lp.,max - 5 X 10^ 1 X IQl" 

ap - 2.5 2.0 

■q = Kp/Ke - 4.8 X 10"* 8 

Up - 1.4 X 10^ 9.2 X 10^ 

Rsrc 1-7 X IQis 5.0 X 10" 4.8 x lO^^ 
B 0.075 80 0.22 

UB 2.2 X lO"'^ 2.5 X 10^ 2 X 10"^ 

Ljet 5 X 10^3 1 X lO^s 5 X 10^^ 

Table 1. The table shows the different parameters used in the modelling of the PKS 2155-304 SED for the three different 
scenarios (SSC, synchrotron proton and mixed scenario) plotted in Figs. 1 to 3. Common values of redshift z = 0.116, 
Doppler factor 5 = 30 and viewing angle 6 — 1° have been used. The normalization parameter K^.p is in units of cm~^, 
and represents the number density of the primary particle distribution at 7 = 1; the size of the emitting ragion Rsrc is in 
cm; the magnetic field in gauss. The energy densities Ue,p,B are given in ergs cm""^; the jet luminosity is in ergs s~^. In 



the evaluation of the jet luminosity the cold proton content of the jet has been included (following Sikora et al. ( 2009 1) 



The authors wish to thank Anita Rcimcr for providing the latest version of the SOPHIA code , and Helene Sol for very useful 
discussions. 
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